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We study a pair of quantum dot exciton qubits interacting with a number of fluctuating charges 
that can induce a Stark shift of both exciton transition energies. We do this by solving the optical 
master equation using a numerical transfer matrix method. We find that the collective influence 
of the charge environment on the dots can be detected by measuring the correlation between the 
photons emitted when each dot is driven independently. Qubits in a common charge environment 
display photon bunching, if both dots are driven on resonance or if the driving laser detunings 
have the same sense for both qubits, and antibunching if the laser detunings have in opposite signs. 
We also show that it is possible to detect several charges fluctuating at different rates using this 
technique. Our findings expand the possibility of measuring qubit dynamics in order to investigate 
the fundamental physics of the environmental noise that causes decoherence. 


I. INTRODUCTION 

Quantum dots (QDs) are semiconductor heterostruc¬ 
tures exhibiting electronic confinement in all three spatial 
dimensions. As such, a QD is zero-dimensional, and its 
eigenstates resemble those of a particle in a box.^^ They 
are ‘artificial atoms’ and many properties typical of a dis¬ 
crete energy level spectra have been observed, for exam¬ 
ple Rabi oscillations .1^1® Qubits may then be represented 
through a variety of different kinds of particles in QDs, 
including electron and hole spin or exciton states.^^^^H^ 
Such QD exciton qubits have large transition dipoles and 
interact strongly with an optical field, and therefore QDs 
make excellent sources of single photons 

A key challenge in assessing the feasibility of any quan¬ 
tum computer realisation is to develop an understand¬ 
ing of decoherence in the system. Though a great deal 
of work has been done on how individual qubits suffer 
decoherence, less is known about how correlated noise 
across multiple qubits arises from different kinds of en¬ 
vironment. Any long-range interaction with an environ¬ 
mental disturbance will cause correlated noise channels 
for relatively closely spaced qubits. One such example 
is that of fluctuating charges in the vicinity of QDs.^^^J^ 
Semiconductors, by their nature, have a Fermi energy in 
the band gap, which can be small enough in doped sam¬ 
ples that the conduction band can be thermally occupied 
when the temperature is relatively low. Vacancies or im¬ 
purities in the crystal structure of a semiconductor lead 
to local alterations to the band structure, and charges 
can become trapped in lower lying states. Depending 
on temperature, such charge traps will randomly switch 
between empty and full, and they may then modelled as 
two state fluctuators. In the vicinity of QD excitons, such 
fluctuators lead to random telegraph noise in the exciton 
energy, and so in t he em itted photon frmuency, through 
the DC Stark shiftsee Figs, [l] and 

In this paper, we show that the common, correlated, 
noise that is generated by charges fluctuating in the vicin¬ 
ity of two QDs can be detected by driving the QDs op¬ 
tically and analysing the subsequently emitted photons. 
Specifically, we determine the cross correlation function 


+ r) of the emitted light, and we will show that 
it can reveal a wealth of information about the nature of 
the charge environment, including how common it is to 
both qubits. In some cases we will find that it is possi¬ 
ble to determine how many charges interact with the QD 
states, and at what rate they fluctuate. 
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FIG. 1: The qubit energy shift for an unoccupied (left) and 
charged (right) trap. The exciton creation energy is denoted 
by uj and the charge-qubit interaction strength is 6. 



FIG. 2: Demonstrative graph of telegraph noise. The 
fluctuating charge switches stochastically between states 1 
(charged) and 0 (uncharged). 

In the following section, we present the model of the 
QDs coupled to both a photonic bath and a common en- 
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vironment of charge fluctuators, and describe the Marko¬ 
vian master equation for the QD-radiation coupling. We 
also present the transfer matrix method for solving the 
time evolution of the density matrix. In section |III| we 
introduce the photon correlation function for two driven 
QDs. In section m we present out predictions before 
summarising in section [V| 


II. MODEL 

We consider two driven QDs, that are modelled as two 
level systems with different energy spacings. These emit 
photons of different frequencies, which can be measured 
by two time-resolving detectors. The QDs are driven by 
lasers with different frequencies to match their respective 
resonances, though each can be slightly detuned from this 
condition. In addition, these uncoupled QDs both inter¬ 
act with a common environment that takes the form of a 
limited number of charge fluctuators that will also be rep¬ 
resented by two level systems; the situation is shown in 
FigU The fluctuators will be treated as classical objects 
with no coherence between the charged and uncharged 
states. The charges affect the qubits via the Coulomb 
interaction, inducing a DC Stark shift. We will neglect 
any dielectric screening effects, which would anyway sim¬ 
ply introduce a renormalisation of the effective distance 
between qubits and fluctuators. 



FIG. 3: Diagram of the experimental set up. The two QDs are 
driven by lasers of different frequencies. The emitted photons 
are captured by the detectors, where the Hanbury Brown- 
Twiss style experiment is performed. 


The Hamiltonian is written as (for full details of the 


Hamiltonian construction see Appendices and : 

i k j 

+ X Ci,k(akO-_,i + a^a+^i) 
k,i 

+ X COs{uJiit) 

i 

+ ( 1 ) 

hj 

The first term represents the individual energies uji of the 
QDs and is the 2 ;-Pauli spin operator for QD i. The 
second term represents the photon energies 6>k of the bath 
of photons with wave vectors k; a\^ and Uk are photon 
creation and annihilation operators The next term rep¬ 
resents the energies of the classical charge traps with 
rji = -\-l for an occupied trap and —1 for an unoccupied 
trap. The next Jaynes-Cummings terms represent the 
interaction of the qubits with the bath of photons with 
coupling strength Q^k- The next term describes the effect 
of two lasers, one coupled to each QD, where uju is the 
laser frequency and Vti is the coupling to QD i. The final 
term is the interaction 5ji between the charge j and QD 
i] Sji = 1 when the charge trap is occupied, and is zero 
otherwise. 

We model our system by using a separate 4x4 density 
operator for each possible charge configuration, each of 
which evolve under the Hamiltonian parameters for that 
particular configuration. By stacking these on top of each 
other, we effectively create a rectangular matrix descrip¬ 
tion of our system, as discussed in Appendix B: Thus al¬ 
together we have an object with 4 x (4n) elements where 
n is the number of possible charge trap states. When a 
charge fluctuates, we simply swap the density operators 
into a new order, reflecting the changed charge configu¬ 
ration. If the charges traps were quantum objects then 
we would require (4n) x (4n) elements to represent it - 
so we save a factor of n in the state description by doing 
this. 

In order to proceed we treat the bath of photons as 
weakly coupled to the QDs and express the dynamics of 
the QDs - fluctuators system in terms of a density matrix, 
leading to the quantum optical master equation,!^ 

+ Y,'nW + i)r[a.,i,pit)] 

i 

+ ^^iNiT[a+,i,p{t)] (2) 

i 

where Hi includes all but the Jaynes-Cummings and pho¬ 
ton energy terms in Eq. following a transformation into 
a frame rotating with the two laser frequencies and a ro- 
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tating wave approximation: 

* j 

+55 + 5Z (^) 

* 

where Ui = uji — uju are the laser detunings. The dissipa- 
ters are given by 

T[L,p{t)] = (Lp{t)L^ - \DLp{t) - \p{t)DL^ . ( 4 ) 

The parameters in Eq. are 7 ^, which is the optical de¬ 
cay rate for the QD i and Ni^ which is the Bose-Einstein 
occupation number of the photon bath taken at the tran¬ 
sition frequencies of QD i. 

Our calculation proceeds using a transfer matrix ap¬ 
proach. This is performed by first writing Eq. in the 
form 

p{t) = Mp{t), ( 5 ) 

where M is the super-operator acting on the density op¬ 
erator p. Since we are solving the system numerically 
and for small time steps, we can write 

p{t + At) = (1 + MAt)p(t). ( 6 ) 

This tells us that the transfer matrix for a system that 
does not fluctuate, but has a constant charge bias can be 
defined as 1 + MAt. We can now introduce the fluctuat¬ 
ing nature of the charge by using a Pauli-x operator, a^c 
(where the ‘c’ indicates that it operates on a charge) to 
flip the charge between its two states (see Appendix [B|) . 
We then divide the time steps into two parts: we assume 
that there is a probability PAt within the time At that 
there is a change in the fluctuator’s state, and a prob¬ 
ability that there will be no change 1 — PAt^ In this 
paper, we will restrict the model to having the same rates 
for hopping onto and away from the trap sites correspon¬ 
dence to a temperature higher than the energy gap of the 
fluctuators. If At is small then higher order terms can be 
neglected and we find that the (1 + MAt) term in Eq. 
is modified to: 

(1 — PAt^]\dAt^PAtdxc ~ 1 T (Af — P-\-P(jxc)At. 

( 7 ) 

Eq. defines the transfer matrix for fluctuating charges 
and with this we can calculate the dynamics of the sys¬ 
tem. 


III. CALCULATING INTENSITY 
CORRELATIONS 

Erom the dynamical simulations of the system density 
matrix, we obtain predictions of emitted photon correla¬ 
tions. To this end we calculate the two-photon intensity 


correlation function 

( 2 )^ {a\{t)al{t + T)a 2 {t + T)ai{t)) 

(a| {t)ai {t + r)a2 {t + r)) 

The operators and at refer to the photon field de¬ 
tected by a detector i G 1,2 (see Eig. [^. Each detec¬ 
tor is responsive only to a range of frequencies around 
the resonant frequency of each QD, and for QDs that 
are sufficiently detuned from one another we expect that 
photons which activate detector i originate from QD i 
only. The QDs in a typical experimental set up would 
be quite closely spaced and so a Hanbury Brown-Twiss 
set upPIl could be used. In this experiment, the photons 
from the QDs are passed through a beam splitter so that 
there are two paths that each lead to a detector, before 
which filters and polarisers can be placed depending on 
what exactly the experiment requires. The detection of 
a photon in the first detector begins a timer, which is 
then stopped by a detection of a photon in the second 
detector.^ The results are collected into a histogram to 
display the number of events as a function of time be¬ 
tween detection events. 

The cross-correlation in Eq.j^is in terms of photon cre¬ 
ation and annihilation operators, but we can relate these 
field operators to our system operators through input- 
output theory.l^ In general an output field is a sum of 
contributions from an input field and from the decay 
of the systems (QDs in our case) that decay optically: 
aout{t) = ain{t) + y^cr_,i(t) + y^cr_, 2 (t). We presume 
a typical setup in which a^n is in the vacuum state, and 
so the output field is then the sum of two well frequency- 
resolved fields, which will be separately detected. We 
can then associated each detector field with a particular 
system operator: 


ai{t) = (9) 

This leads to an expression for g^‘^^(Lt+T) which depends 
only on the QD system operators 

\ = (c ^+,1 + r)a-,2it + T)a-p{t)) 

( 10 ) 

A general way of finding the two time correlation func¬ 
tion from the master equation is to exploit the quantum 
regression theorem.l^ However, we can user a simpler 
method if we assume that experimental measurements 
take place when the system has reached a steady state 
(ss) with density operator pss- In this case is only 
dependent on the delay time r: 

^ (ff+,i(0)o-+,2(r)cr-,2(r)(T-,i(0))ss 

(^+,1*^—,l)ss (<^+, 2 ^^—,2)55 

tr[g+,1(0)0-+,2(r)cr-,2(r)o--,i(0)pss] 

(o^+,lO"—,l)ss(o"+,20'—,2)55 
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Owing to the cyclicity of the trace we can project the 
steady state density matrix into the ground state: 

tr[c7+,2(T)(T_,2(r)(i(l|pss|l)i)|0)ii(0|] 

^ u —w:;—^ 

\cr+pcr_p;ss(cr+^2cr_,2/ss 

where |l)i and | 0 )i represent the two basis state vectors 
for QD 1 . We may now define a new density operator pp, 
which is properly normalised and represents the steady 
state projected from the excited state of QD 1 to its 
ground state: 


i(l|/o..|l)i|0)ii(0| 

,l)ss 


|0)ii(0| 


and write 


(13) 


^ tr[(T+,2(r)g-,2(T)pp] ^ tr[(T-,2pp(T)g+,2] 
(0'+,20- —, 2)55 (cr+,2Cr_ , 2 ) ss 

(14) 

After a long enough period, pp{r) becomes pss and so 
limT-^oo[^^^H^)] = 1 as expected. 


IV. RESULTS 

In this section we present the results across a 

variety of parameters and for a single fluctuator or two 
fluctuators. 


A. One charge fluctuator 

We begin with the case of a single fluctuator. Our 
aim is to assess what kinds of photon cross correlation 
signatures are obtained for different charge fluctuation 
rates and interaction strengths. As a starting point, 
we will assume there is a single charge fluctuator that 
affects each of the two qubits in the same way - i.e. 

= ^12 in Eq. We use QD parameters typical of 
InGaAs structures: Fixed throughout the paper will be 
the spontaneous decay rates (71 = 72 = 1 GHz), and 
the Rabi frequencies (Ui = U 2 = 1 GHz). Other param¬ 
eters are varied for particular sets of results but their 
default values will be: Laser detuning = P 2 = 0, 
charge fluctuation rate is P = 1 MHz. Using values of 
0.8 nm for the permanent dipole and —34 nm^/V for 
the polarizability,!^ for our default choice of charge in¬ 
teraction strength pn = pi 2 = 1 GHz, we would need a 
charge at a distance of 1.32 pm. As the charge is brought 
closer to the qubits it detunes the qubits either towards 
or away from the laser frequency depending on the initial 
size and direction of detuning. 

In Fig[^ we show as a function of pn = pi 2 

between 0 and 10 GHz. This corresponds to effective de¬ 
tunings of 0 and 6 peV, which encompasses the photolu¬ 
minescence range from resonance to effectively zero pho¬ 
ton emission for InGaAs type QDs.!!^ We can see that the 


greater the charge-qubit interaction, the greater the ini¬ 
tial cross-correlation of the detected photons. This initial 
correlation then decays back to the no-correlation value 
of = 1. For all values of the interaction strength 
this decay is on the 1 ps scale, which corresponds to 
the charge fluctuation rate. The explanation for this is 
straightforward: for a larger interaction, then either QD 
is only likely to emit when the charge trap is empty. If 
one QD emits, then since the noise is correlated, the other 
is likely to also emit - at least for over a timescale less 
than the charge fluctuation time. On the other hand, for 
a smaller interaction - less than the Rabi frequency Q 
- then it is also possible for a QD to be excited when 
the charge trap is occupied: an so we expect no cross 
correlations for pn <C 



FIG. 4: cross-correlation for two qubits interact¬ 

ing with a single charge, where the charge-qubit interaction 
strength (Sij) is being varied. 7 ^, the spontaneous decay rate 
of the qubits, is set to 1 GHz and determines the timescale of 
the return to uncorrelated photons. 


We next keep the charge-qubit interaction constant, 
and in Fig. look at how a changing fluctuation rate af¬ 
fects g^^\r). When the fluctuation rate is smaller than 
the photon emission rate, we simply find that g^‘^\r) de¬ 
cays on a timescale similar to that of the charge fluc¬ 
tuation rate. At the fastest fluctuation rates studied, 
however, there is a decrease in the initial value p^^^(O); 
this happens when the charge fluctuation rate exceeds 
the photon emission rate. In this regime, the experiment 
is no longer sensitive to the charge fluctuator. 
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FIG. 5: cross-correlation for two qubits interacting 

with a single charge, where the fluctuation rate (P) of the 
charge is being varied. The values of P are given relative to 
ji (the spontaneous decay rate of the qubits). 

Since the effect of a charge fluctuator is to shift the 
resonance frequency of the qubits away from that of the 
lasers, we can use laser detuning as a further probe of 
the fluctuator correlation dynamics. In Fig. we illus¬ 
trate a typical photoluminescence spectrum of two QDs. 
We know that if the two lasers are resonant with the two 
QDs then we expect to see correlated emitted photons as 
shown in Fig.[^ Imagine instead allowing both the lasers 
to be detuned from the QD resonance (see Fig.[^, by the 
same amount and in the same direction. Depending on 
the direction of the detuning, the effect of the charge will 
be to bring the QDs back into resonance or take them 
further from resonance. In this way, the emitted pho¬ 
tons, regardless of the direction of detuning, will display 
a positive cross correlation. 


Emission Spectrum 

Resonance at Resonance at 

OJl 0)2 



FIG. 6: Schematic drawing of the shifting detuning caused by 
charge interactions when the lasers are detuned in the same 
direction. 

On the other hand if the lasers are detuned in opposite 
directions, i.e. if one is blue shifted and the other red 
shifted with respect to the QD resonant frequencies, then 
charging the trap will shift one QD towards resonance 
and the other further away from it. This situation is 
shown in Fig[^ In this way the QD closer to resonance 
is more likely to emit a photon, while the other is less 
likely. We then expect a negative correlation between 


the emitted photons, for a large enough initial detuning. 


Emission Spectrum 


Resonance at Resonance at 

(jJi 0)2 



FIG. 7: Schematic drawing of the shifting detuning caused by 
charge interactions when the lasers are detuned in opposite 
directions. 


T (^s) 



FIG. 8: cross-correlation for two qubits interacting 

with a single charge, where the laser detuning (yi) is being 
varied. On the left of the zero detuning line, detuning in¬ 
creases and is equal in magnitude and sign for both dots; on 
the right it increases with equal magnitude but opposite signs. 


In Fig we show the cross-correlation for both the 
detuning scenarios just described. As expected, we see 
positive correlation for same sense detuned, and nega¬ 
tive for opposite sense - expect around z/ = 0, when the 
charge fluctuation moves both QDs away from resonance 
and we recover a positive cross-correlation. 


B. More than one charge 

Let us now introduce a second charge fluctuator into 
the model, and establish whether it is possible to dis¬ 
tinguish multiple from single fluctuators by using cross 
correlation measurements. 
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FIG. 9: cross-correlation for two qubits interacting 

with two charges, where the fluctuation rate (Pi) of the 
charges are being varied. The values of P are set equal to 
each other. 
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FIG. 10: cross-correlation for two qubits interacting 

with two charges, where the fluctuation rate (Pi) of one charge 
is being varied and the other held at 10“^ GHz. 


Fig.j^shows g^‘^\T) as a function of (equal) fluctuation 
rates for the two charges. We can see immediately that 
g^‘^^ (0) is larger in this case than for the single fluctuator 
case shown in Fig. This is expected since the two 
fluctuators working together increase the total possible 
detuning of each QD; as we have seen in Fig.j^this results 
in a higher initial cross-correlation. As the fluctuation 
rate increases, the decay of g^^\r) happens at shorter 
times, similar to the single charge case. 

It is unlikely, however, that the two charge traps will be 
fluctuating at exactly the same rate, so let us now look 
into how different rates of fluctuation affect the cross¬ 
correlation. In Fig. 10, we show the impact on g^‘^\r) of 


altering one fluctuation rate while keeping the other fixed 
at 10“^ GHz. Comparing this figure with that for a single 
fluctuator shown in Fig. we find that in certain cases 
it is possible to see a clear qualitative difference between 
the results for a single and two fluctuators. This is easier 
to discern by taking cuts through the plots for particular 
fluctuation rates, and using a log scale for the time; for 
the single and two fluctuator cases these are shown in 
Fig, prj^a) and Fig. [T^b). If the rates for the two fluctu¬ 


ators are significantly different then two plateaux can be 
seen in the curves, with two decay rates corresponding 
to two different fluctuation rates; this effect washes out 
once the faster fluctuation rate approaches that of the 
QD optical decay rate. 


a) 

g®(T) 



FIG. 11: Log time plots for some representative frequencies 
of the data shown in Fig. shown in a) and Fig. 10 shown in 
b). 


As would be expected, when the fluctuation rates of 
the charges are very similar, no deviation from the single 
fluctuation curves can be distinguished: They do not ex¬ 
hibit the plateaued structure that can be seen in FigpT] 
b). 

Finally, we look at varying the laser detunings for two 
qubits and two charges. In Fig we show the cross cor¬ 
relation function for the same detuning parameters as in 
Fig-i in the case where the two fluctuation rates are not 
equal (1 GHz and 1 MHz). We also display various cuts 
through this 3D plot, for different values of the detuning, 
in Fig. There is a clear contrast here with the sur¬ 
face shown in Fig. the cuts in Fig. [T^ show that the 
plateaued structure characteristic of two different rates 
survives as detuning is varied. At negative detuning the 
a negative correlation is observed, but with a clear long 
time plateau. The oscillatory behaviour as detuning in¬ 
creases is a consequence of the increasing effective Rabi 
frequency, which is no longer fully damped by the 1 GHz 
optical decay processes. 
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FIG. 12: g^‘^\T) cross-correlation for two qubits interacting 
with two charges, where the laser detuning {vi) is being varied. 
The two charges fluctuate with different rates of 1 MHz and 
1 GHz. 





FIG. 13: Guts through Fig.[^ for five values of the detuning 
parameters. 


V. SUMMARY 

We have described a measurement which can reveal in¬ 
formation about the long range correlations in the deco¬ 


hering environment of a controlled quantum system. We 
have shown that for an environment consisting of one or 
two fluctuating charge traps that the time dependence 
and initial value of the cross-correlation generate a sig¬ 
nature of a correlated environment and in some cases can 
give an indication of how many charge traps are present. 

Figs. [8| and [T^ display the variation of the experimen¬ 
tal signature as a function of the detuning. Since detun¬ 
ing can be varied in a single experiment, these are key 
predictions that provide the most experimentally acces¬ 
sible signature of correlation, and clear differences are 
observed for the single and two charge trap cases. The 
switching from a positive to a negative correlation as a 
function of detuning serves as a distinctive signature of 
a common environment; the double plateaued structure 
is a sign of two distinct fluctuation rates. 

As the numbers of charges increase beyond two, sev¬ 
eral changes are expected to the plots. The initial 

correlation of the emitted photons will be higher, due to 
the larger number of charges giving an effective greater 
charge interaction strength in line with Fig. Addition¬ 
ally, if the charges have widely varying fluctuation rates, 
then each charge would be seen as a separate plateau. 

Experiments which detect the effects we have described 
would not rely on any fast detectors or particularly fast 
optics. They are relatively simple measurements that 
could be performed on a suitable sample immediately, 
and would provide a unique probe of the unexplored col¬ 
lective effects of open system environments. 
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a 4 X 8matrix, can be treated as two 4x4 matrices, each 
corresponding to one of the states of the charge. 

With this in mind we define the system Hamiltonian 
as 



The first column vector describes the charge state; uji is 
energy required to excite uncoupled QD i from ground 
to excited state. ^ is the energy difference between the 
charge being at the trap site and the charge being else¬ 
where, which we assume to be at an infinite distance. 
Each operator is now of dimension 4x8; during coher¬ 
ent evolutions the top and bottom square matrices of 
this object are treated individually and act on the top 
and bottom square matrices of the corresponding density 
operator respectively, following the usual rules of quan¬ 
tum mechanics. Operator eigenstates corresponding to 
the different charge configurations can be found for each 
square component matrix. 


Appendix A: Qubits and Charge Hamiltonian 

Here we present the Hamiltonian construction in the 
case of two QDs and a single charge; more complex sit¬ 
uations follow similarly. Each QD i {i G {1,2}) has 
a ground state |0)i and excited state as \l)i, so that 
their uncoupled Hamiltonian is denoted in a spin lan¬ 
guage as (Jz^i = |l)i(l| — |0)i(0|- The charge fiuctuator 
two level systems are classical and so cannot have coher¬ 
ences. While a quantum two level system is represented 
as a two by two density matrix, a classical object has 
no coherences, so the off-diagonal matrix elements are 
unnecessary. A two-vector of classical populations (rjj) 
then describes the state of fiuctuator j, with the first el¬ 
ement giving the population in the ground (no-charge) 
state, and the second that in the excited (charged) state. 
If we use the usual tensor product formulation of quan¬ 
tum mechanics of more than one subsystem, our overall 
description of our qubit-fluctuator state must then be a 
rectangular matrix, corresponding to 2^ stacked square 
matrices, one for each of the classical states of N fiuc¬ 
tuator s. Between charge fluctuation events each square 
density operator acts independently of the others and so 
each can be treated individually. Eor example, the case 
of two qubits and a single charge, which is described by 


Appendix B: Interaction Hamiltonian and 
Fluctuations 

The qubit-bath interaction Hamiltonian is defined as: 

Hq^B = Cl,k(<^k^-,l +Q^k^+,l) +C2,k(<^k^-,2 +<^k^+,2), 

k 

(Bl) 

where 

<7+,i = ^ j ^ (g)l 2 (g)cr+, 

0’+,2 = ^ J ^ (g)cr+ (g)l2, 

(T_,i = (^ j ^ ® I2 ® (T-, 

(T-^2 = (^|^®(T+®l2 . (B2) 

The interaction here describes the couplings between the 
QDs and the photonic bath, with coupling strengths 
or C 2 ,k, where Uk and a^. are the annihilation and creation 
operators for the photon bath mode with wave vector k. 
The QD raising and lowering operators in the interaction 
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Hamiltonian are denoted by the cr+ and cr_ Pauli matri¬ 
ces, which are constructed as |1)(0| and |0)(1| respec¬ 
tively. The effect of this part of the interaction with the 
system is taken into account by deriving a Born-Markov 
optical master equation. 

We also account for a coupling between the qubits and 
the (classical) lasers. The qubit-laser Hamiltonian is 


Hq^l = ftlCFx,! COs{uJiit) Q2(Tx,2 COs(cJ/ 2 t), (B3) 

where Qi is the Rabi frequency of the ith QD and uju is 
the laser frequency of laser i which is assumed to drive 
only the ith QD with which it is closely resonant. Simi¬ 
larly to previous definitions: 







( 8 ) 12 ® (Jx, 


l 2 - 


(B4) 


The final interaction is between the qubits and charge: 

IT _ ^11 I ^12 /-OrN 

Hq^c — 

where 

Ain = (^ g ^ 0 I 2 

IJ'12 = (^0^00-^012. (B6) 

There are two terms, which show the interaction of the 
charge (labelled 1 here since we are considering only a 
single charge) with both of the QDs. In general the inter¬ 
action is where the subscript denotes an interaction 
of charge j with QD i. The term describes a process in 
which the stationary charge at some well defined distance 
introduces a Coulomb potential that gives rise to a Stark 
shift of the excitonic states, thus creating a TLS with a 
larger energy spacing. 

The qubit-laser Hamiltonian is time dependent, but 
we can remove that time dependence by moving into a 
rotating frame and performing the rotating wave approx¬ 
imation (RWA). This is done by applying a unitary trans¬ 
formation to all part of the Hamiltonian that involve the 


dots and the charges, i.e., all parts except Hq^s- From 
the time dependent Schrodinger equation, we get: 

- ihdtiu-^ I -4’)) = -ih{dtU-^) I -ip) + u-^H I -ip) 

= HU-^ I p}) 

Therefore, the effective Hamiltonian becomes, after mov¬ 
ing to the rotating frame 


= UHU-^ + U{ihdtU-^). (B7) 

The operator U~^ is 

/ 0 0 0 \ 

0 g 5 *t(wi 2 -Wil) Q Q 

0 0 Q 

\ 0 0 0 j 

(B8) 

Notice that this is a 4 x 4 matrix. As stated previously, 
we can consider the two charge states as independent and 
as such, the 4x8 rectangular matrix can be thought of as 
two 4x4 matrices. This unitary transformation operates 
on each of those two parts individually. After performing 
the rotating wave approximation, which removes rapidly 
oscillating terms, i.e. those with twice the frequency of 
ujii and uji 2 ^ and letting uji — uju = z/^, we get: 


H' = I 


+ ^2<^z,2 + + ^2<^x,2) 

1 ^ 


(B9) 


This is now our final system Hamiltonian; in the paper 
we show how this is taken together with Hq^B to find 
a quantum optical master equation for the qubit-charge 
system. 

Additionally, we need to define an operator in which 
a single charge fiuctuates from occupied to unoccupied 
or vice versa^ an operation executed by the operator cFxc- 
In our rectangular density operator notation for a sin¬ 
gle charge this simply corresponds to swapping over the 
upper and lower square matrices. For more complex sit¬ 
uations of more charge fiuctuators, a similar reordering 
of the now multiple square matrices achieves the desired 
effect. 



